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Abstract 

The Higgs sector of the [/(l)-extended supersymmetric model is studied with 
great detail. We calculate the masses of the Higgs bosons at the one-loop level. We 



■ also calculate at the one-loop level the gluon-involving processes for the productions 

43 . and decays of the scalar Higgs bosons of the model at the energy of the CERN 

Large Hadron Collider (LHC), where the radiative corrections due to the loops of 
CN ■ top, bottom, and exotic quarks and their scalar partners are taken into account. 

We find that the exotic quark and exotic scalar quarks in the model may manifest 
£N) \ themselves at the LHC, since the production of the heaviest scalar Higgs boson via 

gluon fusion processes is mediated virtually by the loops of exotic quark and exotic 
scalar quarks, for a reasonable parameter set of the model. 

O 

O- 1. Introduction 



The search for the Higgs boson is one of the prominent goals of the future high-energy 
collider experiments. The massive fermions and gauge bosons in the standard model 
(SM) acquire their masses through the mediation of a single scalar Higgs boson. The 
mass of this SM Higgs boson receives quadratic divergences at the one-loop level due to 
the radiative corrections from other SM particle loops. The simplest way to solve this 
problem is to introduce the supersymmetry (SUSY) to the SM. The SUSY partners to 
the SM particles help to eliminate the quadratic divergences to the Higgs boson mass, 
by the cancellation between each SM particle loop and its superparticle loop, at the one- 
loop level. However, if the SUSY is a broken symmetry in nature, the cancelation is not 
complete and thus the Higgs boson mass receives non-zero quantum corrections at the 
one- loop level [1-4]. 

The most economical SUSY model is the minimal supersymmetric standard model 
(MSSM), where two Higgs doublets are introduced in order to cancel the gauge anomaly. 
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They give masses independently to down-type fermions and up-type fermions. Since 
studies of the MSSM has revealed a few theoretical problems such as the \x problem [5], 
various alternative SUSY models have been introduced. 

The U (l)-extended supersymmetric standard model (USSM) is one of the non-minimal 
supersymmetric models which does not suffer the MSSM fi problem. The USSM is char- 
acterized by the extra U(l)' gauge symmetry at the TeV scale [6-23]. Thus, there is an 
extra neutral gauge boson with a mass of a TeV scale. Also, the USSM has a larger Higgs 
sector, which consists of a neutral Higgs singlet besides two Higgs doublets. In the USSM, 
the quantity corresponding to the \x parameter of the MSSM is dynamically generated by 
the vacuum expectation value of the Higgs singlet. 

Another characteristic of the USSM is that each SM generation is extended by an extra 
pair of SU{2) singlet quarks, D l L and D l R , with electric charges —1/3 and +1/3, respec- 
tively. The presence of these exotic quarks is interesting from the phenomenological point 
of view. In particular, scenarios with a light exotic quark in the USSM have extensively 
been investigated [18-22]. A recent study have shown that the exotic quarks may con- 
tribute by large amount to the CP mixing between the heaviest scalar and pseudoscalar 
Higgs bosons [18]. Since the exotic quarks couple directly to the neutral Higgs singlet, 
the heaviest Higgs boson mass at the one-loop level might be affected by the exotic quark 
loops. Also, the exotic quark loops contribute to the gluon fusion process of the Higgs 
production at the Large Hadron Collider (LHC). These effects of the exotic quarks have 
been studied elsewhere [22]. While some approximations have been made in Ref. [22], 
the calculations in this article are carried out in a rigorous way. 

In this article, we continue to study the Higgs potential in the USSM at the one-loop 
level. In particular, we are interested in the Higgs productions and decays at the LHC. 
We take into account the radiative corrections due to the quarks of the third generation, 
exotic quarks, and their superpartners at the leading order. We calculate the production 
cross sections of the scalar Higgs bosons via the gluon fusion process. It is found that 
the production of the heaviest scalar Higgs boson might be significantly enhanced for a 
parameter region by the exotic quark and scalar quarks contributions. 

We also study the dominant decay modes of the scalar Higgs bosons. Among them, 
we calculate the decay processes of the scalar Higgs bosons into gluon pairs, which are 
mediated by the loops of the quarks, exotic quarks, and their superpartners. The gluonic 
decay mode might be important for the Higgs search of the USSM at the LHC, since it 
may exhibit the role of the exotic quark and exotic scalar quarks. In particular, we find 
that the decay of the heaviest scalar Higgs boson into a pair of gluons depends virtually 
only on the loops of the exotic quark and exotic scalar quarks. 

The article is organized as follows: In the next section, the Higgs potential of the 
USSM is described. Then, we calculate the masses of the pseudoscalar and the three 
neutral scalar Higgs bosons at the one-loop level. These one-loop masses are taken as 
input for the Higgs production calculations in Section 4 and the Higgs decay calculations 
in Section 5. Numerical analysis is presented in Section 6, and conclusions in the last 
section. 
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2. Higgs Potential 



There are two Higgs doublet superfields Hi and H 2 , and a Higgs singlet superfield S in 
the USSM. The gauge symmetry of the USSM is 577(2) x U(l) x U{1)'. For the Higgs 
superpotential, we consider only the quarks of the third generation, and one family of the 
exotic quarks, Dl and D R . Thus, the superpotential of the USSM may be written as 

W « h t Q T eH 2 f R - h h Q T eUxb c R + h k SD L D R + \SH\eH 2 , (1) 

where h t , hb, and h k are the Yukawa couplings of the top, bottom, and exotic quarks, 
respectively, Q is the left-handed quark doublet superfield of the third generation, t c R and 
b c R are respectively the charge conjugate superfields of the right-handed top quark and 
the right-handed bottom quark, A is the dimensionless coupling constant, and e is a 2 x 2 
antisymmetric matrix defined as €i 2 = 1. 

The Higgs sector of the USSM consists of two Higgs doublets, Hi = (Hf,H~) and 
H 2 = (H + , H®), and one Higgs singlet S. The tree-level Higgs potential, V°, is obtained by 
collecting the F-terms and the U-terms of the above superpotential, and adding physically 
allowed soft terms. Thus, it may be written as 

V° = V F + V D + U s , (2) 

with 

Vf = \\?[{\H^ + \H^)\S\ 2 +\HleH^\ , 

V D = 9 4(HtaH 1 + HlaH 2 f+ 9 4(\Hi\ 2 -\H 2 \ 2 r 

o o 

+ 9 f(Qi\Hi\ 2 + Q 2 \H 2 \" + Q 3 \S\ 2 )\ 
U s = m\\Hi\ 2 + m\\H 2 \ 2 + ml\S\ 2 - [\A x H? e H 2 S + H.c] , (3) 

where g 2 , gi, and g[ are respectively the gauge coupling coefficients of SU{2), U(l), 
and U{1)', a are Pauli matrices, A\ is the trilinear soft SUSY breaking parameter with 
mass dimension, rrtj (i = 1,2,3) are soft SUSY breaking masses, and Qi, Q 2 , and Qs are 
respectively the effective U(l)' hypercharges of Hi, H 2 , and S. We do not consider in 
this article the CP violation in the USSM, which may take place in the mixing between 
the scalar and pseudoscalar Higgs bosons. Thus, the parameters of the Higgs potential 
as well as the vacuum expectation values of the neutral Higgs fields are taken to be real 
in following procedure. 

The effective U(l)' hypercharges of the Higgs fields satisfy the identity relation of 
J2i=i Qi — in order to ensure the U{1)' gauge invariance. The soft SUSY breaking masses 
rrii may be eliminated by using the three minimum conditions for the Higgs potential 
that define the vacuum. After electroweak symmetry breaking, the neutral components 
of Higgs fields will develop vacuum expectation values as v\ = (Hf), v 2 = (H®), and 

s = (S). We introduce a free parameter, tan/3 = v 2 /v\. The value of v = \Jv 2 + v 2 is 
fixed by electroweak data as 175 GeV. In terms of the vacuum expectation values, the 
masses of top, bottom, and exotic quarks are respectively given as m t = h t v 2 , m h = hf,Vi, 
and mfc = hkS after electroweak symmetry breaking. 
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The tree-level 2x2 mass matrices for the scalar partners to top, bottom, and exotic 
quarks in the USSM may be written as 



M q n Mf 2 



M q = —11 —12 (4) 



V Mil M% 2 

where q = t,b,k. For simplicity, we rewrite Mfj as Mfj, Mfj, and Mjj, respectively, for 
the scalar partners to top, bottom, and exotic quarks. The matrix elements are obtained 
from the tree-level Higgs potential in the on-shell Lagrangian after the elimination of the 
auxiliary fields as follows: 

Ml, = mJ + / i M 2 +(f - !) (I^l 2 " |tf 2 T) 
+ ^QQiQ^l 2 + Q 2 \H° 2 \ 2 + Q 3 \S\ 2 ) , 
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,2 



M\ 2 = m 2 u + h 2 AH 2 2 \ 2 + 9 -±{\Hl\ 2 -\H»\ 2 ) 

+ ^QuiQ^l 2 + Q 2 \H° 2 \ 2 + Q 3 |5| 2 ) , 
M{ 2 = h t (\H°*S* - A t H 2 ) , 

+ 9 -Iq q (Q x \hI\ 2 + Q 2 \H° 2 \ 2 + Q 3 \S\ 2 ) , 



1 



,2 



M 2 = m 2 D + K 2 h \Hl\ 2 -^{\Ht\ 2 -\Hl\ 2 ) 

o 

+ ^Q^Q^ 2 + Q 2 \H° 2 \ 2 + Q 3 \S\ 2 ) , 
M\ 2 = h b (\H$*S* - A h H») , 

= m 2 K + h 2 k \S\ 2 + 9 ^\ 2 -\Hl\ 2 ) 

+ ^Q^Q^l 2 + Q 2 \H° 2 \ 2 + Q 3 \S\ 2 ) , 
M| 2 = ml + ^-^d^T-l^l 2 ) 



'2 



(3 



+ ^q^q^ 2 + g 2 |i/ 2 T + g 3 |5| 2 ) , 

M\ 2 = h k (\H°*H° 2 * - A k S) , (5) 

where mg, mjy, mp, and are the soft SUSY breaking masses for the scalar quarks, 
A t , Ab, A k are the trilinear soft SUSY breaking parameters for them, and Qq = —Qi/3, 
Qu = (Qi - 2<52)/3, Qd = (Qi + 2Q 2 )/3, Q K = Q2, and = Q\. These relations for 
the effective U(l)' hypercharges are dictated by the cancellation of anomalies. Note that 
the contributions from Vd are included in the above mass matrices for the scalar quarks. 
In Ref. [22], they have been neglected. Thus, if we set g\ = g 2 = g[ = 0, the above 
expressions will reduce to the corresponding results in Ref. [22]. 
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Now, the squared masses of the scalar partners to top, bottom, and exotic quarks 
are respectively given by the eigenvalues of (q = t, b, k), by diagonalizing them. We 
obtain the squared masses of the scalar quarks as 



in 



ml i 

0l,0 2 



m Q + m2u + m 2 t + ^f cos 2(3 + G[v 2 cos 2 + G\v 2 sin 2 f3 + G\Q,s< 



™q + m 2 D 
2 

m 2 K + m 2 ^ 



m 2 h -^ cos 2/3 + G\v 2 cos 2 p + G b 2 v 2 sin 2 (3 + G b 3 s 



b n 2 



+ m'l + G\v 2 cos" (3 + G h 2 v A sin" (3 + G*s" T 



fe„.2 -2 



(6) 



where the contributions of the scalar quark mixing are given by 



X h 



Xh 



m 2 Q -m1j 



+ 



2m 2 y 5m 2 



cos 2/3 + B\v 2 cos 2 /3 + Bl sin" (3 + B\s 



t „2 



+ m 2 (Xs cot (3 — A t 



rrin — m 



D 



m K — m 



K 



K 3 


12 


- A t ) 2 


• 


' ' m 2 z 


m w\ 


J2~ 


3 ) 


-A b ) 2 


J 


(m\ 


m 2 ^ 


3 


3 , 



>b„,2 



?6„,2 



?b„2 



?fc„2 



+ m 2 , ( At> 2 sin /3 cos f3/s — A k 



(7) 



with 



G< = 






G\ = 




+ Q D )Qi 


G\ = 


9 f(Q K 


+ Qr)Q 1 


B\ = 


'2 


— Qu)Qi 


B\ = 


g f(Q Q 


— Qd)Qi 


B\ = 


9 4(Qk 


— Qk)Qi 



3. Higgs Masses at the One-Loop Level 

At the one- loop level, quarks and their scalar partners contribute to the radiative cor- 
rections through loops. The contributions from the loops of top quark and top scalar 
quarks are most dominant for a wide region in the parameter space. Also, the radiative 
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corrections due to the loops of bottom quark and bottom scalar quarks might be phe- 
nomenologically significant, especially for very large tan/3, at the low energy scale. The 
contributions from the loops of the exotic quark and its scalar partners in the USSM are 
not well studied. We have suggested elsewhere that their contributions are worthwhile 
studying from a phenomenological point of view because the exotic quark couples directly 
to the Higgs singlet S in the USSM [22]. Thus, in this article, we consider the loops of 
top, bottom, exotic quark and their scalar partners in our calculations at the one-loop 
level. 

The Higgs potential at the one- loop level, V 1 , is given as 



V 1 = V° + Vl 



where the radiative corrections, Vi jE ff, is obtained by effective potential method [24]. 
Explicitly, it is given by 



v - V n * M * 



64tt 2 



M] 3 



(9) 



where A is the renormalization scale in the modified minimal subtraction scheme, and 
n q = — 12 is the quark degree of freedom and = 6 (i = 1, 2) are the degrees of freedom 
for scalar quarks, determined from their color, charge, and spin factors. 

After the electroweak symmetry breaking, ten real degrees of freedom from the Higgs 
sector of the USSM are reduced to six physical Higgs particles, namely, a pair of charged 
Higgs bosons, one pseudoscalar Higgs boson, and three scalar Higgs bosons. The squared 
mass of the neutral pseudoscalar Higgs boson at the one-loop level may be written as 

m A — m A° + m A* + m A b + m A k > (10) 
where m^o is the tree-level mass, given as 

2 2\vA x 

m A° = ■ , 
sin 2a 

and mA*, ^A h i an d m A k account respectively for the radiative corrections from the top, 
bottom, and the exotic quark sectors: 

2 3mf\A t , 2 2 \ 

m At = z o~7, t\ m i , m f ) , 

A 87^ sin 2 /3 sin 2a v tv t2 > ' 

2 3m b \A b ft 2 2 \ 

77T. 4b = 7. / 1711 , TTlr , 

A 8vr 2 t;cos 2 (3 sin 2a V 61 ' h *> ' 

3m 2 k XA k v 2 2 

^ f(m~ , ml) , (11) 

87r 2 s 2 sin2a' /v k ^ ' v 1 



m 2 Ak 



where a is a mixing angle between the electroweak scale and the U (1)' symmetry breaking 
scale, defined as 

v 

tana = — sin 2(3 , (12) 
2tS 
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and the dimensionless function /(mf , m|) is defined as 



[mj — mf) [ A 2 A 2 



(13) 



The three neutral scalar Higgs bosons, Si (i = 1, 2, 3), of the USSM are given by the 
eigenvectors of a 3 x 3 symmetric mass matrix, and their squared masses, rris i (i = 1, 2, 3), 
are given by the corresponding eigenvalues. These neutral scalar Higgs bosons are sorted 
such that m Sl < m S2 < m S2 . 

At the one-loop level, the mass matrix for the three scalar Higgs bosons, may be 
decomposed as 



where is the matrix elements at the tree level, obtained from V°, given explicitly as 



and M-j, M\j , and Afg are respectively the radiative contributions from the top quark 
sector, the bottom quark sector, and the exotic quark sector. They are obtained from 

v i,eff as 




m 2 z cos 2 + 2g' 2 Q\v 2 cos 2 /3 + m 2 A o sin 2 j3 cos 2 a , 

2 -2010 /2-o2 2 • 2 a i 2 2 a 2 

m z sm p + 2p 1 Q 2 v sin p + m A o cos p cos a , 

<-> /2/S2 2 i 2-2 

V3 S + m A° sm a > 

g'iQiQ2 v2 sin 2/3 + (A 2 t> 2 — m z /2) sin 2/3 — m^o cos /3 sin /3 cos 2 a , 
2g' 2 QiQ 3 vs cos /3 + 2A 2 t> s cos /3 — m 2 A o sin /3 cos a sin a , 
2g l2 Q 2 Qzvs sin /3 + 2\ 2 vs sin /3 — m^o cos /3 cos a sin a , 



(14) 




327r 2 -y 2 
3 



32vr 2 t; 2 




) 



(15) 



where q = t,b,k and 



2 

(2G> 2 + ^)cos/3 , 





(2G^ + ^)sin/3, 



2G\vs , 
2GX cos/3 , 
2G^ 2 sin/3 , 




(16) 
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-i 2m 2 \sA t , , 9 4mL 5mi. ^ A „ 
= \ R * + (2Blv 2 + — 2- cos/3A? , 
sin p 3 6 

= - 2m ^f * - (-2i^ 2 + ^ - ^1) sin/3A 
2 sin/3 v 2 3 6 ; ' 

w , = 2niXvA, +2B , nsAlt 



tan (3 

2mlA b A h o 2m?* 



W? = -^^ + (25 1 V-^ + ^)cosM, 
cos p 3 6 

cos p 3 6 

W 3 ft = 2m 2 Av tan (3A b + 2B b 3 vs A 9 b , 

W? = M^^ + (25y + ^i-^i)cos/3Af , 

W* = 2^A^ S cos^A fc _ ( _ 2BV + 2m| _ 2m| ) ^ 
s 3 3 

,2 



2 



W* = _ 2 - 4 W Afc + 2j gfc i;gA g ; (17) 



A t = As cot /3 -A t , 

4 5 

Af = mg-my + ( - m 2 ^ - -m%) cos 2/3 + 2E^ 2 cos 2 (3 + 2fi2^ 2 sin 2 /3 + 2.63s 2 , 

o 
Ab = As tan j3 — Ab , 

1 2 

A£ = -m 2 D + (-m? z - -m 2 w ) cos 2(3 + 2B b v 2 cos 2 /3 + 25^ 2 sin 2 /3 + 25^s 2 , 

6 3 

At — Awtano — yl^, (18 
A' t = m|. - mjj + -(m| - mfy) cos 2/9 + 2Bft> 2 cos 2 /3 + 2Bjt> 2 sin 2 /3 + 2Bf s 2 , 

CS, = ,^sm 2 /!cos 2 a-^(2BX + ^-^) 2 /K,<), 

Z4 = r^,cos^cos^-^(-2By + ^-5|i) 2 /K, m|) 

3m^ f m t\ 
~ 4n 2 v 2 sin 2 /3 ° g V A 2 "/ ' 

£>*3 = m^sin 2 ^- i^-/( m ?, m|) , 

. 9 „ 3 sin 2/3 / , 9 4mL 5m 2 7 \ 

D\ 2 = - m 2 ,co S ^ s in/?c OS 2 a + ^3^- 2Sy + ^£--^ 



x(-2^ 2 + ^-^)/K,<) 



*1> "^2' 
o l,OS /3 2 2 



3m 2 A 2 s cos/3 

L) 13 = - m At sin £ cos a sin a — = — r^-rj{m- t m t 

8ir z v sin p 



3j5oScos/3 / „ + 4m 2 ,, SmiX „. , 
' 25^ 2 + ^ ^ f\m\ , m 2 



^7T 2 f \ 3 6 
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l% 2 v \ 2 3 6 / 1 2 



3m ' bg ( <\ 

4tt 2 v 2 cos 2 /3 V A V ' 
3 sin 



n () 2 2 a 2 3 shl 2 /3 / . 2 2m^ . 2 2 \ 

L> 22 = m A6 cos /3cos a - -2.5ft; ^ + ^ /( m 6i> m &J 



^33 = m V in2a -l^-/K' m U 



3 sin 2/3 /__>, 2 2mL mj 
L>?2 = - m 2 4 i cos 8 sin 5 cos 2 a H ^ 2RV 

327T 2 V 2 1 



-2B 2 V-^ + ^) /(m|, m|) , 



-D 13 = — m A b sin 8 cos a sin a 

3B^scos8 ( nrth 2 SmL m 2 z \ „. 2 2 . 

' 2B\v 2 - + f(ml , ml) , 



87r 2 w V 1 3 6 / 6l ' 62 
n b 2 o 3m 2 A 2 stan/3 2 2 
^23 = - "i^b cos ^ cos a sin a — = —}{m~ b mj 



35oSsin/3 / 6 2 2mL m|\ , 2 



57r 2 v I 2 3 6 / v 61 ' b2 ' 



D k 


= m 2 A h sin 2 8 


D k 
^22 


= m\ k cos 2 8 


n k 
-^33 


= m 2 A k sin 2 a 


u \2 


= — m A h cos 



1% _ 2m| N 
I 3 
3sin 2 /3 / n , 2 2m| 2m w 



3cos 2 /3 ( nr>k 9 2mi 2?r% , 

— £ 25> j H - f(m( , mi , 

16n 2 v 2 I 1 3 3 / kl fc2 



167rV \ 2 3 3 J J v fc2 

2 \ 
it 1 



(^) 2 ,,2 2 v 3mt , /m 2 



3^/ 2 2^_2^' 
32vr 2 t; 2 \ 1 3 3 

3ml\ 2 v 2 sin 2/3 „. 9 9 . 

-f{ml, m 2 



16tt 2 s 2 ' v fey ' 

3m?A 2 s cos /3 

8n 2 vsm 2 8 J K kl k2 ' 



r^k 2 ■ n • * JI " J k awo M p / 2 

iJiq = — m 4fc sm p cos a sin a o-^rt\ m T - mr 



3Bkscos8 f nT , k o , 2m z 2m w \ 
^ — - 2B k v 2 H — f(m( , m( ) , 

-D 23 = — m^fe cos /3 cos a sin a 
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+ + ^ - ^) /(mi, <) . (21) 



The above analytic expressions are obtained by employing some mathematical formu- 
las. We remark that, in a previous study, we have neglected for the sake of simplicity the 
.D-term contributions of the the scalar quarks for the neutral Higgs boson masses [22]. 
One may obtain the results of Ref. [22] by taking g\ = g 2 = g[ = in the above expres- 
sions. Thus, the above expressions for the scalar Higgs boson masses at the one-loop level 
are more rigorous and complete, and calculated for the first time in the USSM. 



4. Higgs Productions in pp Collisions 



It is well noticed that the gluon fusion process, mediated by the top quark triangular 
loops, is the most dominant production channel for the SM Higgs boson in pp collisions 
at the energy of the LHC. The gluon fusion process is also important for the productions 
of the scalar Higgs bosons of the USSM at the LHC, where the loops of top quark and 
top scalar quarks, by far, would play the most dominant role, for a large space of the 
relevant parameters of the USSM. For very large tan/3, the loops of bottom quark and 
bottom scalar quarks would also yield considerable contributions. Further, the loops of 
the exotic quark and exotic scalar quarks may also contribute as much significantly as the 
loops of top quark sector or the bottom quark sector. Thus, in our analysis, we consider 
all of them, namely, the contributions of top, bottom, and the exotic quark sectors to the 
productions of the USSM scalar Higgs bosons via the gluon fusion process. 

Let us denote &i(gg —> Si) as the parton- level cross section for the gluon-gluon anni- 
hilation into Si (i = 1,2,3) in the USSM. It is obtained as [25,26,22] 



Oi(s) = &i{gg -> St 



® 2 s (™z) 
256tt 



m- 



m 



St 



A S D \ 2 ) 



(22) 



where s is the square of the cm. energy of two gluons, a s (mz) is the strong coupling 
constant evaluated at the electroweak scale, and A s represents the contributions from 
top and bottom quarks and their scalar partners while A S D accounts for the contributions 
from the exotic quark sector. They are given in terms of form factors as 



A' 



A s 



E 

q=t,b 



G Siqq A s q (T q ) 



3=1,2 



G Sl kkA s k (T k ) + J2 G s t k,k 

3=1,2 



V 



J k J m 



2 



(23) 



where A q (r q ) is the quark form factor (q = t, b, k), Gs m is the coupling coefficient of Si 
to a pair of quarks (q = t, b, k), A~ ijq) is the form factor of the scalar quarks (q = t, b, k), 
and Gsiqjq-j are the coupling coefficients of 5^ to the pairs of scalar quarks (j = 1,2 and 
q = i,b, k). 
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The quark form factor and the form factor for the scalar quarks are given respectively 

as 



A S q {r q ) = rjl + (1 - r q )f{r q )} , Af (r qj ) = r qj [r qj f(r qj ) 
where the scaled variables are defined as 



1] 



Am 2 



Am 2 - 



m 



Si 



m 



(24) 



(25) 



Si 



and the function / is defined as 



f(r) 



arcsin 2 (l/ v /r) 



Ml 



l +yl+T 



Z7T 



T > 1 
T < 1 



(26) 



Now, let us calculate the coupling coefficients, Gs m and Gs iqjqj - It is convenient to 
normalize the coupling coefficients of Si to a pair of quarks using the corresponding SM 
coupling coefficient. Then, the normalized coupling coefficients for our case are obtained 

as 



G f 



cos (3 



G 



SM 



sin/3 



(27) 



where are the elements of the 3x3 orthogonal matrix which diagonalizes the mass 
matrix for the scalar Higgs bosons. The coupling coefficient of Si to a pair of exotic 
quarks, normalized with the exotic quark mass m^, is given as 



G 



-0 M . 

s 



The coupling coefficients of Si to a pair of scalar quarks are defined as 

3 



G 



Siqj<lk 



i=i 



(2« 



(29) 



where r 5 ^~ (i = 1, 2, 3 and q = t, b, k) are the 2x2 matrix for the coupling coefficients 
of Si to a pair of the left- and right-handed components of weak eigenstates of scalar 
quarks, Z7* is the 2x2 orthogonal matrix which diagonalizes the mass matrix for the 
scalar quarks, and the subscript indices are /3,7 = L, R, and j, k = 1,2. Note that 
are also present in the above expression. 

The explicit expressions for T Si ^ are as follows: 



iSitt 




-,9i 





'2 

-ttQuQi 



v 2 sin (3 



11 



pS3&6 



+ u sin /3 



6 2 8 

v 



2 
3'' 



~<7l + ^0^2 



m t I Acos/3 \ 
■usin/3 \ Acos/3 y 



'2 

-ttQuQ?, 



2m h -A h 



v 2 cosP \ -A b 2m b 



+ v cos j3 



( 



\ 



2 2 '2 
91 92 , £1 ^ 



6 



m fe / As 
v 2 cos/9 V As 

+ v sin j3 



\ 



91,92,91 
6 2 











4 + €qdQ2 



Asin/3 



ucos/3 \ A sin ,5 



+ s 



V 







'2 

^QdQ 3 



m k I Asin/3 
s I A sin /9 



+ v cos f3 



1 g * + 9 '*n n 



V 







9l < 9i SS 



QkQ 



1 J 



m k I Acos/3 
s \ A cos f3 



+ u sin (3 



( 



\ 



2 '2 
^ + -£-QkQ2 
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The 2x2 diagonalizing matrix for the scalar quarks, U^, is given as 



U q 



cos 9„ 



sin 9 a 



sin 9a cos On 



(31) 



where the mixing angle Q„ is defined as 



tan 9„ 



Mf 2 



(32) 



with being the mass of the lighter scalar quark and M|- being the elements of the 
mass matrix M q for the scalar quarks (q = t, b, k). 

Now, we are ready to calculate the production cross section of the scalar Higgs bosons 
at the LHC. In order to obtain (Ji{pp — » Si), the desired cross section at the laboratory 
level, we should fold the cross section at the parton level with the gluon distribution 
functions as [25,26,22] 



r 1 r 1 fg 

ai(E,m Si )= dr dx—g(x,m s .)g(r/x,m 

J Tg J T X 



<a As 



tE 2 



(33) 



where g(x,m s .) is the gluon distribution function at the factorization scale m s ., with x 
being the fraction of the momentum of the proton carried by the gluon, g(jsl 'x,m s .) is 
the gluon distribution functions for the other participating gluon, and r s = m s ./E 2 is the 
Drell-Yann variable, with E being the cm. energy of the colliding protons at the LHC. 
For the gluon distribution functions, the CTEQ6L is used in our calculation [28]. 



5. Higgs Decays 

The scalar Higgs bosons would decay, among others, into pairs of fermions. In the SM, the 
pairs of bottom quarks would be the most prominent decay modes for the SM Higgs boson 
if its mass is below 135 GeV. One may sort out the branching ratios for the SM Higgs 
boson with an intermediate mass of around 110 GeV as BR(Hsm bb) > BR(Hsm 
tt) > BR(Hsm — > gg). Thus, for such an intermediate SM Higgs boson, decaying into 
a pair of gluons via the triangular loop of the top quark would be more difficult than 
decaying into a pair of bottom quarks or a pair of tau leptons. On the other hand, if the 
mass of the SM Higgs boson turns out to be larger than 135 GeV, the Higgs decay into 
a pair of gauge bosons would be more dominant than any other decay modes. 

The situation in the USSM, we expect, would be different from the SM case, since 
there are exotic quark and exotic scalar quarks in the USSM. In particular, they would 
certainly affect the decays of the scalar Higgs bosons. As one can see in the superpotential 
of the USSM, the heaviest scalar Higgs boson couples directly to the exotic quarks, while 
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the other scalar Higgs bosons couple indirectly to the exotic quarks through the mixings 
among them. Thus, in order to study the effects of the exotic quarks in the decays of 
the scalar Higgs bosons, we may as well pay more attentions to the heaviest scalar Higgs 
boson than the other ones. Let us study the decay modes of the scalar Higgs bosons of 
the USSM into pairs of fermions, gluons, and gauge bosons, as they are considered to be 
the most important decay modes. Among them, the decay mode into gluon pairs should 
be studied in more detail, since it may reveal the effects of the exotic quarks more clearly 
than other decay modes. 

The total decay width of the Si (i — 1,2, 3) may be assumed as 

r(Si) = T(Si -> bb) + T{Si -> rr) + T{Si -> /i/i) + T(Si -> cc) 

+ r(S, -> ss) + r(3i -> gg) + r(S, -»• WW) + r(S, -> ZZ) , (34) 

where each decay mode can be understood without difficulty Note that the heaviest 
scalar Higgs boson, S3, has additional decay modes, since it can decay into the lighter 
ones, that is, S3 — > Si Si, S 3 — > S 2 S 2 , or S 3 — > SiS 2 . 

Now, let us describe each partial decays. The partial decay width of Sj into a pair of 
fermions is given as 



C f mjm s , 

167TV 2 



Am 



Gsiff » ( 35 ) 



where Cf is the color factor of fermions, m/ is the mass of the fermion, and Gsjf is the 
coupling coefficient normalized with respect to the corresponding SM coupling coefficient. 
This fermionic partial decay width applies to bottom, charm, and strange quarks, and 
tau and muon leptons. 

The partial decay width of S, into a pair of gauge bosons is given as 



+ 12-f G% vv , (36) 



m" q \ m"q m% 



where Cz = 1 and Cw = 2, m v is the mass of the gauge boson, and GsiW is t ne coupling 
coefficient of Si to the gauge boson pairs, normalized by the corresponding SM coupling 
coefficient (V = Z,W). 

The normalized coupling coefficient of Si to a pair of Z bosons is given as [21] 

G Sl zz = cos PO^I + sin f30 2l C 2 2 + ^^C 3 2 , (37) 

^HZZ 

where Grzz = gi^zl cos#iy is the corresponding SM coupling coefficient, with 9w being 
the weak mixing angle, and C, (z = 1,2, 3) are dimensionless parameters defined as 

C\ = cos 2 9 W + sin 2 9 W cos (ft — — — - - cos 9 W sin 9 W sin <\> , 
C2 = cos 2 9w + sin 2 6w cos <j) + ® x ^ 2 cos 9w sin 6w sin <fi , 

C3 = g'lQz sin 9 W sin <p , (38) 
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with <p being the mixing angle between Z and the extra neutral gauge boson in the USSM. 
Notice that <fi in Cj induces the interference effect of the extra neutral gauge boson. 
Also, the normalized coupling coefficient of Si to a pair of W bosons is given as 

GsiWW = cos (50 u + sin (50 2i . (39) 

Note that it does not depend on <p, since there is no extra charged gauge bosons in the 
USSM. 

The partial decay width of Si into a pair of gluons is given as 

m^ g9 ) = a ' { ™£f* (\AS + AS D \*), (40) 

where mediating loops involve top, bottom, and exotic quarks and their scalar partners. 

Now, as we have mentioned before, the heaviest one of the three scalar Higgs bosons 
may decay into the lighter members. The partial decay width for Si — > SjSk is calculated 

as 



T(S l SjS k ) -- 



(m s - m s J 2 \ / (m s .+m Sk y 



(41) 



where Cs is the symmetry factor and Gs j s k Si is the cubic Higgs coupling coefficient. We 
have Cs = 1 for j ^ k and Cs = 1/2 for j = k, and 

G SjSjSi = X 2 s{20 ll O u 3l + 2 lj 3l + 20 2l 2l 3j + 2 2j 3l ) 
- XA x (O u 2j 3j + O^OMj + O l3 2j 3t ) 
+ v\ 2 cos $[0^0% + Oy + 20 l3 (0 2j 2l + 3j 3i )} 
+ v\ 2 sin (3[0 2l (0 2 + Ol) + 20 2j {0 lj O u + 3] 3i )} 



+ -Ml + <? 2 ) cos /3[30^0 l4 - O h O^ - 20 l3 2j 2l )) 



'lj ^ ^3j 

(n 2 4- cos 

r 
4 



+ V -{g 2 + ^ 2 2 ) sin f3 [30 2 2j 2l - 2l 2 Vj - 20,,0 2i O H )] , 

+ ^ 2 g 3 s[20 1 ,Oi i o 3i gi + o 2 ,o 3 iQi + 20 2 ,o 2i o 3j Q 2 
+ o|o 3i g 2 + 30 2 3j o 3i Q 3 } 

+ g' 2 v cos PQi^OlpuQi + O u O%Q 2 + 20 lj 2j 2i Q 2 
+ O u 2 3j Q 3 + 20 l3 33 3l Q 3 ] 

+ g' 1 2 v S ml3Q 2 [20 lj 1 i0 2j Q 1 + O^Qi + 30*^4 
+ 2% 0%Q 3 + 20 2 ,0 3 ,0 3i Q 3 ] , 

G Sl s 2 s 3 = - ^XA x (0 13 22 31 + 12 23 31 + 13 21 32 

+ Ou0 23 32 + o 12 o 21 o 33 + o n o 22 o 33 ) 

+ X 2 s(0 12 13 31 + 22 23 31 + 11 13 32 

+ 2 i0 23 32 + O n 12 33 + 2l 22 33 ) 

+ vX 2 cos f3[0 13 2 i0 22 + 12 2 i0 23 + On0 22 23 

+ Oi 3 3 i0 32 + Oi 2 3 i0 33 + O n 32 33 ] 

+ vX 2 sin f3[0 12 13 21 + O u 13 22 + O u 12 23 

+ o 23 o 31 o 32 + 22 31 33 + 21 32 33 ] 
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+ 4(^1 + 9l) cos (3[30 u 12 13 - 13 21 22 

— 0\ 2 2 \0 23 — 0\\0 22 23 \ 

+ \{g\ + al) sm (3[30 21 22 23 - 12 13 21 

- O u 13 22 - O u O l2 23 \ 

+ ^ 2 g3s[0 12 13 3 iQi + 11 13 32 Q 1 + 11 12 33 Q 1 

+ 22 23 31 Q 2 + 21 23 32 Q 2 + 21 22 33 Q 2 + 30 31 32 33 Q 3 ] 

+ g*v cos PQi [30 n 12 13 Qi + 13 21 22 Q 2 + 12 21 23 Q 2 

+ O n 22 23 Q 2 + O13O31O32Q3 + 12 31 33 Q 3 + O11O32O33Q3] 

+ ^?uBin^Q 2 [0 12 1 30 21 Q 1 + 11 1 30 22 Q 1 + 11 12 23 Q 1 

+ 30 21 22 23 Q 2 + 23 31 32 Q 3 + 22 31 33 Q 3 

+ 21 32 03 3 g3] • (42) 



6. Numerical Analysis 

Now, we are ready for carrying out numerical analysis. The scheme for setting the relevant 
parameters is similar to our previous analysis [18-22]. First of all, we note that there are 
strong experimental constraints in the USSM on the mass of the extra gauge boson, mz>, 
and on the size of the mixing angle \<p\ between Z and Z'. In this article, we will use 
rriz' > 700 GeV and |0| < 3 x 10 -3 . The vacuum expectation value of the Higgs singlet 
is set as s = 700 GeV in order to ensure a light Z 1 scenario with a mass below 1 TeV. 
The effective U(l)' hypercharges of the Higgs doublets and Higgs singlet can be redefined 
as Qi = g[Qi because Qi appear always together with g\. Then, we fix the redefined 
effective hypercharges as Q\ = —1, Q 2 = —0.1, and Q 3 = 1.1. The Tevatron Run I data 
have set the lower bound on the mass of an long-lived quark with electric charges ±1/3 
as 119 GeV [29]. The search for long-lived charged massive particles at Tevatron Run I 
put a more stringent experimental lower bound of 180 GeV at the 95 % confidence level 
[30]. Thus, we take the mass of the exotic quark as 400 GeV without contradicting the 
Tevatron constraint. Also, the mass of the pseudoscalar Higgs boson at the one-loop level, 
rriA, will be used as an input parameter. 

The present experimental lower bound on the mass of the SM Higgs boson set by LEP2 
data is about 114.5 GeV [31]. Since the scalar Higgs bosons in the USSM are mixtures 
of the real components of the three Higgs fields, the possibility of discovering them is 
affected by how much the real component of the Higgs singlet is mixed in them. Further, 
the decay channel for the scalar Higgs bosons in the USSM is not restricted to a pair 
of bottom quarks. Thus, the search for a USSM Higgs boson would depend on not only 
its mass but also its coupling coefficients and other factors. Therefore, in our numerical 
analysis, the lower bound on the mass of the lightest scalar Higgs boson in the USSM may 
not be restricted by the LEP2 constraint, although it does not imply that our numerical 
analysis may phenomenologically contradict the experimental constraint of LEP2. Our 
analysis have been performed by assuming that the scalar Higgs bosons decay exclusively 
to a pair of bottom quarks. 

It is well known that the top scalar quark plays a dominant role in the Higgs processes 
involving gluons, namely, both the productions of a scalar Higgs boson via the gluon 
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fusion process at the LHC and its decays into a gluon pair, in the parameter region where 
the top scalar quark mass is below 400 GeV [25,26]. The result depends crucially on 
the mass of the top scalar quark. In the USSM, as we concentrate on the role of the 
exotic quark and the exotic scalar quarks, the masses of the exotic scalar quarks should 
be carefully studied. Thus, it is of importance to investigate the dependence of the exotic 
scalar quark masses on the trilinear soft SUSY breaking parameter, Ak, and on the mass 
of the exotic quark, m^. 

We show the result for the exotic scalar quark masses as functions of Ak in Fig. 1, 
for tan/3 = 3, A = 0.387, mg = mu = m D = m K = = 300 GeV, A t = A b = 440 
GeV, and rrik = 400 GeV and = 150 GeV. The mass of the lighter exotic scalar quark 
is calculated to be between 176 GeV and 428 GeV while the mass of the heavier one is 
between 493 GeV and 628 GeV. Their masses should be degenerate if Ak vanishes, but 
actually they are not exactly degenerate in mass since the non-zero A invokes a mixing 
between them. 

Note that the masses of the other scalar quarks as well as the scalar Higgs masses do 
not depend on Ak- The masses of the other scalar quarks are calculated as m~ tl = 229 
GeV, mi 2 = 435 GeV, = 252 GeV, and m^ 2 = 339 GeV. The masses for the three 
scalar Higgs bosons at the one- loop level are calculated as ms 1 = 117 GeV, ms 2 = 156 
GeV, and ms 3 = 783 GeV, which are also shown in Fig. 1. 

In order to study the effects of the exotic quark sector in the production of the scalar 
Higgs bosons in the USSM via the gluon fusion process at the LHC, it would be convenient 
to calculate the production cross sections with or without the contributions of the exotic 
quark sector. Let us denote the production cross section of Si as cr h (Si) where only 
the exotic contributions are taken into account, and the production cross section of Si 
as a t,b (Si) where they are not included. These cross sections are obtained technically 
by taking either A s = or A S D = in the formula for the production cross section. 
We also introduce a (Si) where the loops of top, bottom, and exotic quarks and their 
scalar partners are all included. It should be noted that a(Si) is not exactly the sum of 
& t,b (Si) and a k (Si), since there are interference terms in + Af,| 2 in the formula for the 
production cross section. However, the interference between them are negligible, for the 
parameter values we consider, such that a (Si) is almost equal to a t,b (Si) + a k (Si). The 
results are shown in Figs. 2a and 2b, where the values of the relevant parameters are the 
same as in Fig. 1. 

In Fig. 2a, a t,b (Si) and a k (Si) are plotted as functions of Ak. Note that a t,b (Si) 
do not change against Ak because both the scalar Higgs boson masses and the relevant 
couplings are independent of Ak, whereas the dependency on Ak of a k (Si) is very clear. 
Note also that a t ' b (S 3 ) is negligibly smaller than either a t,b (S2) or a t,b (Si). This is mainly 
because the mass of S3 is much heavier than the other two scalar Higgs bosons. The trend 
is consistent with the SM prediction that the production cross section decreases as the 
Higgs mass increases. It is also worthwhile noticing that a k (S 3 ) is about 10 3 times larger 
than a t,b (S 3 ). This fact indicates that the loops of exotic quark and exotic scalar quarks 
virtually dominate the S3 productions through the gluon fusion process in pp collisions. 

In Fig. 2b, a(Si) (i — 1, 2, 3) are plotted as functions of Ak- Comparing with Fig. 2a, 
one may note that cr(Si) = a t,b (Sx), cr^) = cx*' 6 ^), but 17(83) — c k (S3), in particular 
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if Ak is not so much large as 400 GeV. Therefore, we may expect that the productions 
of the heaviest scalar Higgs bosons of the USSM at the LHC would definitely exhibit the 
contributions of exotic quark and exotic scalar quarks in the USSM. 

Now, let us study decay processes. We calculate the partial decay width of Sj (i = 
1, 2, 3) into a pair of gluons. Here, too, it would be convenient to separate the contributions 
of exotic quark and exotic scalar quark loops from the contributions of top and bottom 
quark and scalar quark loops. Let us denote the partial decay width of Si into a gluon 
pair as T k (Si) where only the exotic contributions are taken into account, and the partial 
decay width of Si as r* ,6 (Sj) where they are not included. We also introduce r(Sj) where 
the loops of top, bottom, and exotic quarks and their scalar partners are all included. 

We show r t,fe (5'j) and T k (Si) in Fig. 3a. They are plotted as functions of A^. The 
other parameter values are the same as in Fig. 1. Like a t,b (Si), one can see that T t,b (Si) 
are all constant against A k . Note that r t,6 (5i) > r*> 6 (S 2 ) > r*' b (S 3 ), but the difference is 
not significant. On the other hand, there is a great hierarchy among r fc (S;). Note that 
r fc (S , 3) is larger than either r fc (Si) or r fc (S 2 ) by several orders of magnitude. 

In Fig. 3b, r(Sj) (i = 1,2,3) are plotted as functions of A^. One may notice that 
r(Si) 7^ r t,6 (S;) + T k (Si). This is due to the interferences between them. Clearly, by 
comparing Fig. 3a and Fig. 3b, the discrepancy between T(Si) and r(S 2 ) is more visible 
than the discrepancy between r*' 6 (Si) and r t,6 (S 2 ). We expect that the interference 
between r*' 6 (S 2 ) and r fc (S 2 ) is rather destructive such that r(S 2 ) < r t,fe (S 2 ). Nevertheless, 
for the heaviest scalar Higgs boson, it is evident that practically r(Ss) = r fc (S , 3), since 

r fc (S3)»r'> fe (s 3 ). 

Next, we turn to the branching ratios of Si (i = 1,2,3). The results are shown in 
Figs. 4-6. The values of the relevant parameters are the same as in Fig. 1. In Fig. 4, 
we show the branching ratios of Si, namely, BR(Si —> bb), BR(Si — > cc), BR(Si — > rr), 
and BR(Si — > gg), as functions of A^. The other branching ratios are very small to be 
included in the figure. Note that the branching ratios for pairs of ordinary quarks are 
flat against Ak, whereas BR(Si — >■ gg) increases very slightly as A k increases, since the 
mass of the lighter exotic scalar quark decreases down to about the top quark mass as A^ 
increases. Comparing Fig. 4 with the following figures, one can see that BR(Si — > bb) is 
much larger than BR(S 2 — > bb), or BR(S3 — > bb) that is too small to be shown in Fig. 6. 

In Fig. 5, we show the branching ratios of S 2 , namely, BR(S 2 —> bb), BR(S 2 — > 
cc), BR(S 2 — > rr), and BR(S 2 — > gg), as functions of Ak- Here, one may notice that 
BR(S 2 -» cc) > BR(S 2 -»■ bb), and BR(S 2 -> gg) > BR(S 2 -> rr), while the corresponding 
inequalities are reversed in case of Si, as shown in Fig. 4. Also, notice that BR(S 2 — > gg) 
decreases mildly as A^ increases by the contribution of the exotic quark and exotic scalar 
quarks for the full range of A k . 

Since S3 is much heavier than Si or S 2 , its decay channels are richer than the lighter 
scalar Higgs bosons. In particular, it may decay into a pair of exotic scalar quarks or into a 
pair of exotic quarks, depending on the masses of participating particles. According to our 
analysis, the masses are obtained asms 3 = 780 GeV, = 400 GeV, and 200 < m k < 450 
GeV. Thus, S3 may decay into a pair of exotic quarks off shell, and into a pair of the 
lighter exotic scalar quarks both on and off shell. 

Let us consider the on-shell decays of S3 into a pair of the lighter exotic scalar quarks 
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for Ak > 150 GeV. Its partial decay width is given by 
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where G\g is given by Eq. (29). This formula is used to calculate the branching 
ratio of S3 into a pair of k\. In Fig. 6, we show BR(S*3 — > kik\), the branching ratio 
for S3 decays into a pair of the lighter exotic scalar quarks, as function of Ak, as well 
as other dominant branching ratios of S3, namely, BR(S 3 — > Si Si), BR(S 3 — > S 2 S 2 ), 
BR(S 3 -> SiS 2 ), BR(S 3 ->■ tt), BR(S 3 -> WW), BR(S 3 -» ZZ), and BR(S 3 ->■ gg). 

It is clear in Fig. 6 that S3 would decay mostly into a pair of lighter scalar Higgs 
bosons: BR(S3 — > Si Si) and BR(S3 — > S 2 S 2 ) are almost equal and dominant over other 
decay channels. Nevertheless, BR(S3 — > fcifci) increases sharply at Ak ~ 150 GeV, and 
becomes quite compatible to them. Note that BR(S3 — > k\ki) reaches its maximum at 
Ak = 270 GeV. One may also notice that the branching ratio for a gluon pair is far less 
than 1 %, and fluctuates visibly due to the Ak dependence of the exotic contributions. It 
is comparable to BR(S3 — > WW) or BR(S3 — > ZZ). However, the S3 decays into a gluon 
pair is almost solely mediated by the loops of the exotic quark and exotic scalar quarks. 

In Fig. 6, one may notice that S3 may decay quite dominantly into a pair of Si, a 
pair of S 2 , a pair of the lighter exotic scalar quarks, a pair of top quarks, or into Si 
and S 2 , as well as a pair of gluons. The branching ratio of S3 into a pair of gluons is 
very small as 0.001 at best. We can obtain the number of raw Higgs events by using the 
integrated luminosity of the detectors such as CMS [32] or ATLAS [33] . For the integrated 
luminosity of 700 fb _1 which is expected to be reached in about ten years at the LHC 
(CMS + ATLAS), about 140 raw Higgs events for S 3 would be produced at the LHC [34]. 
Thus, the decay mode of S3 into a pair of gluons would not be relevant for any practical 
purpose, since the background for such events is large at the LHC. Consequently, other 
decay modes would be more worthwhile studying, or the super LHC (SLHC), an upgrade 
of the LHC, is necessary to study the gluon decays of the heaviest quark in our model. 
Note that the target of the integrated luminosity for SLHC is about 3000 fb~ x during the 
period of its operation [34]. 

The mass of the SM Higgs boson is determined by the vacuum expectation value of 
the electroweak symmetry breaking, in which itself plays a crucial role. Thus, it is natural 
to expect that the SM Higgs boson has a mass of the scale of the electroweak symmetry 
breaking. By requiring that there should be no Landau pole up to the cut off scale, the 
renormalization group equation for the mass of the SM Higgs boson predicts that it is 
smaller than 190 GeV, if the cut off scale is the Planck mass scale. 

If the SM Higgs boson is discovered with a mass of the electroweak symmetry breaking 
scale and if S3 of the present model is about 780 GeV, there would be no serious difficulty 
to distinguish between them. On the other hand, if the SM Higgs boson is discovered with 
a mass comparable to the mass of S3 of the present model, their decay modes may be a 
useful tool to distinguish between them. For the parameter value we consider, the mass 
of S3 is about 780 GeV, and its decay modes are shown in Fig. 6. Let us assume that 
the mass of the SM Higgs boson is also about 780 GeV. With this mass, the SM Higgs 
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boson would decay 60 % in W boson pairs, 30 % in Z boson pairs, and 7 % in top quark 
pairs. On the other hand, as Fig. 6 shows, S3 would decay 45 % into Si Si and 45 % into 
SiS 2 for < A k (GeV) < 130, or 30 % each for S1S1, SiS 2 , and hh for A k ~ 250 GeV. 
Further, as Figs. 4 and 5 show, Si and S2 would respectively decay into a pair of bottom 
quarks by 60 % and 19 %. Therefore, using the decay modes of bottom quark pairs as 
background, assuming 700 fb _1 for the integrated luminosity of the LHC, we would have 
29862 events of S3 — > bb whereas 2362 events for the same decay modes for the SM Higgs 
bosons. In this case, the signal significance is more than 614<r, well above the discovery 
limit. 

Since S3 is much heavier than Si or S2, it may decay into a pair of exotic scalar quarks 
or into a pair of exotic quarks. For the parameter values we take up to now, the mass of 
S3 is 780 GeV, the mass of the exotic quarks is 400 GeV, and the mass of the lighter exotic 
scalar quarks is between 180 and 450 GeV. Thus, for this choice of parameter values, S3 
may decay into a pair of the lighter exotic scalar quarks both on and off shell. On the 
other hand, it is kinematically forbidden for S3 to decay into a pair of exotic quarks on 
shell. 

In order to study the on-shell decays of S 3 into a pair of exotic quarks, we need to 
explore other regions of the parameter space. We would like to study the parameter region 
where m k = 400 GeV and ms 3 is larger than 800 GeV. We find that a representative set 
of the parameter values in this region is tan /3 = 3, A = 0.5, vtlq = my = my> = irtx = 
m R = 500 GeV, A t = A b = 440 GeV, m k = 400 GeV, m A = 900 GeV, and < A k < 900 
GeV. These parameter values yield 176 < m ki (GeV) < 585, 635 < (GeV) < 845, 
m Sl ~ 114.5 GeV, m$ 2 ~ 750.6 GeV, and m s . A ~ 928.7 GeV, satisfying our requirement. 
Fig. 7 shows the masses of these particles as functions of A k . Hence, with this choice of 
parameter values, S3 may decay into a pair of exotic quarks on shell. 

We calculate the production cross sections for the scalar Higgs bosons via gluon fusion 
process at the LHC. The radiative corrections include the contributions from loops of 
top quarks, bottom quarks, exotic quarks, and their scalar superpartners. The results 
are shown in Fig. 8. The parameter values are the same as Fig. 7. We then calculate 
the branching ratios of the scalar Higgs bosons. The results for the branching ratios are 
shown in Figs. 9, 10, and 11, for the same parameter values as Fig. 7. 

These figures may be compared with Figs. 4, 5, and 6, respectively. Consider first 
Fig. 9, which shows the branching ratios of Si. The behavior of Si in Fig. 9 is not so 
much different from its behavior in Fig. 4, mainly because its mass is roughly the same 
in both figures: ms 1 ~ 114.5 GeV in Fig. 9 and ms 1 ~ 117 GeV in Fig. 4. 

By contrast, the behavior of the branching ratios of S2 in Fig.10 is strikingly different 
from that in Fig. 5, since m$ 2 ~ 750 GeV in Fig. 10 and m$ 2 ~ 156 GeV in Fig. 5. In 
particular, the branching ratio for a pair of lighter exotic scalar quarks, BR (S2 — > &1&1), 
can be seen in Fig. 10, whereas it is absent in Fig. 5, as it is forbidden kinematically 
there. The coupling coefficient of S2 to a pair of kik\ depends very severely on A k . Thus, 
the branching ratio for a pair of lighter exotic scalar quarks varies wildly as A k changes. 
Note that S2 cannot decay into a pair of exotic quarks on shell, since its mass is smaller 
than 2m k , in both figures. 
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The heaviest scalar Higgs boson, 53, on the other hand, may decay into a pair of 
exotic quarks on shell, as shown in Fig. 11. One may notice that the branching ratio of 
S3 — > kk is considerably larger than the branching ratios for gauge boson pair or gluon 
pair. Therefore, for some parameter regions where 5*3 is allowed kinematically to decay 
into a pair of exotic quarks, the decays of the heaviest scalar quark into exotic quarks may 
become important to examine our model. Also, the branching ratio for a pair of exotic 
scalar quarks in Fig. 11 varies more wildly than in Fig. 6. 

7. Conclusions 

In this article, we study the Higgs sector of the USSM at the one-loop level. The radiative 
corrections from the loops of top, bottom and exotic quarks and their superpartners 
are calculated for the scalar and pseudoscalar Higgs boson masses, where we include 
the D-term contribution. For the numerical analysis, we take tan/3 = 3, A = 0.387, 
rriQ = mjj = m D = m K = = 300 GeV, A t = A b = 440 GeV, and the mass of the 
exotic quark as = 400 GeV and the mass of the pseudoscalar Higgs boson as tua = 150 
GeV. We assume that CP symmetry is conserved in the Higgs sector of the USSM. We 
find that the exotic scalar quarks become nearly degenerate in mass for small A^. 

For the USSM Higgs phenomenology at the LHC, we study the production processes 
and decay modes of the scalar Higgs bosons that involve gluons. As the scalar Higgs 
bosons may be produced in pp collisions via the gluon fusion process, we calculate the 
production cross sections for the gluon fusion process, considering the loop contributions 
from top, bottom, and exotic quarks and their superpartners. We find that, for the 
heaviest scalar Higgs boson of the USSM, the gluon fusion process is mediated virtually 
only by the loops of exotic quark and exotic scalar quarks. So is the partial decay width 
of the heaviest scalar Higgs boson into a pair of gluons. In other words, practically, the 
heaviest scalar Higgs boson couples to a pair of gluons only through the loops of exotic 
quark and exotic scalar quarks. Thus, the contributions from the loops of top and bottom 
quarks and their scalar partners are negligible in the gluon fusion process of the heaviest 
scalar Higgs boson production as well as in its partial decay width into a gluon pair. 
However, the partial decay width of the heaviest scalar Higgs boson into a pair of gluons 
is less than 1 %, just comparable to its partial decay width into a pair of gauge bosons. 

We should note that our results are obtained for a particular set of parameter values, 
allowing A^ to vary between and 450 GeV. Thus, it may be hard to say that this set 
of parameter values represent the whole parameter space of the USSM. We may have 
to search a much wider region of the parameter space of the USSM in order to obtain 
more predictive results for the contributions of the exotic quark and exotic scalar quarks 
in the Higgs phenomenology. However, qualitatively, our results suggest that the exotic 
quark and exotic scalar quarks play observable roles in the productions and decays of the 
scalar Higgs bosons of the USSM at the LHC. In order to examine the USSM at the LHC, 
in particular to examine the effects of the exotic quark sector in the Higgs physics, the 
heaviest scalar Higgs boson may be one of the best windows, as our calculations clearly 
suggests. 
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FIG. 1: The scalar Higgs boson masses ms i (i = 1,2,3) and the exotic scalar quark 
masses m k . (j = 1,2) are plotted as functions of A k , for tan/3 = 3, A = 0.387, ttlq = 
rau = m D = m K = = 300 GeV, A t = A b = 440 GeV, the mass of the pseudoscalar 
Higgs boson as = 150 GeV, and the mass of the exotic quark as m k = 400 GeV. Note 
that rris 3 > ms 2 > ms 1 are independent of A k , whereas the gap between the masses of 
the two exotic scalar quarks widens as A k increases. 
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FIG. 2a: The production cross sections of the scalar Higgs bosons via the gluon fusion 
process at the LHC are plotted as functions of A k . a t,b (Si) (i = 1,2,3) are the cross 
sections where only the loops of top quark, bottom quark, top scalar quarks, and bottom 



scalar quarks are taken into account, and a (Si 



1, 2, 3) are the cross sections where 



only the loops of the exotic quark and exotic scalar quarks are taken into account. The 
parameter values are the same as Fig. 1. Note that a t,b (Si) (i = 1,2,3) are independent 
of A k . 
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FIG. 2b: The production cross sections of the scalar Higgs bosons via the gluon fusion 
process at the LHC, a(Si) (i = 1, 2, 3), are plotted as functions of Ak- They are calculated 
by taking into account together both the loops of top quark, bottom quark, top scalar 
quarks, bottom scalar quarks and the loops of the exotic quark and exotic scalar quarks. 
The parameter values are the same as Fig. 1. Notice that cr^S^) is virtually equal to 
a k (S 3 ) in Fig. 2a. 
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FIG. 3a: The partial widths of the scalar Higgs boson decays into a pair of gluons are 
plotted as functions of A^. T t,b (Si) (i = 1,2,3) are the partial decay widths where only 
the loops of top quark, bottom quark, top scalar quarks, and bottom scalar quarks are 
taken into account, and T k (Si) are the decay widths where only the loops of the exotic 
quark and exotic scalar quarks are taken into account. The parameter values are the same 
as Fig. 1. Note that r*' 6 ^) (i = 1,2,3) are independent of A k . Note also that T k (S 3 ) is 
much larger than any of T t,b (Si) (i = 1,2,3) whereas r /c (S'i) or T k (S2) are smaller than 
them. 
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FIG. 3b: The partial widths of the scalar Higgs boson decays into a pair of gluons, T(Si) 
(i = 1,2,3) are plotted as functions of A^. They are calculated by taking into account 
together both the loops of top quark, bottom quark, top scalar quarks, bottom scalar 
quarks, and the loops of the exotic quark and exotic scalar quarks are taken into account. 
The parameter values are the same as Fig. 1. Note that is practically identical to 

r fc (s 3 ). 
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FIG. 4: The branching ratios of the lightest scalar Higgs boson, Si, are plotted as 
functions of A^. The parameter values are the same as Fig. 1. Note that BR(Si — > bb), 
BR(S'i — > cc), and BR(Si —> tt) are independent of Ak, whereas BR(Si — > gg) slightly 
increases as Ak increases. 
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FIG. 5: The branching ratios of the middle heavy scalar Higgs boson, £2, are plotted as 
functions of A k . The parameter values are the same as Fig. 1. Note that BR(S < 2 — > bb), 
BR(S 2 — > cc), and BR(S' 2 — > rr) are independent of A k , whereas BR(S , 2 — > gg) slightly 
decreases as A k increases. Note also that BR(S < 2 — > gg) is larger than BR(S , 2 — > tt). 
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FIG. 6: The branching ratios of the heaviest scalar Higgs boson, S3, are plotted as 
functions of A k . The parameter values are the same as Fig. 1. Note that BR(S < 3 — > tt), 
BR(S 3 -> WW), BR(S 3 ->■ ZZ), as well as BR(S 3 -> S 1 S 2 ), BR(S 3 -> Si^), and 
BR(S 3 —7- S , 2 S , 2 ), are independent of A k , whereas BR(S , 3 — > kiki) and BR(S , 3 — > gg) 
depend significantly on A/,. Note also that BR(S' 3 — > gg) is larger than or quite compatible 
to BR(S 3 -> WW) and BR(S 3 -> ZZ). 
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FIG. 7: The same as Fig. 1, except for a different set of parameter values: tan f3 = 3, 
A = 0.5, rriQ = mjj = = = rug = 500 GeV, A t = = 440 GeV, the mass of 
the pseudoscalar Higgs boson as = 900 GeV, and the mass of the exotic quark as 
mfc = 400 GeV. Note that ms 3 is larger than twice the mass of the exotic quark. 
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FIG. 8: The same as Fig. 2b, except that the parameter values are those of Fig. 7 
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FIG. 9: The same as Fig. 4, except that the parameter values are those of Fig. 7 
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FIG. 10: The same as Fig. 5, except that the parameter values are those of Fig. 7 
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FIG. 11: The same as Fig. 6, except that the parameter values are those of Fig. 7. Here, 
note that the branching ratio for a pair of exotic quark pairs, BR(S < 3 — > kk) is present. 
It is the second largest branching ratio, larger than BR(S , 3 — > WW), BR(S , 3 — > ZZ), or 
BR(S 3 ^gg). 
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